On the Chaotic Flux Dynamics in a Long Josephson 

Junction 



Z. C. Feng and Y. Charles Li 



Abstract. Flux dynamics in an annular long Josephson junction is studied. 
Three main topics are covered. The first is chaotic flux dynamics and its pre- 
diction via Melnikov integrals. It turns out that DC current bias cannot induce 
chaotic flux dynamics, while AC current bias can. The existence of a common 
root to the Melnikov integrals is a necessary condition for the existence of 
chaotic flux dynamics. The second topic is on the components of the global 
attractor and the bifurcation in the perturbation parameter measuring the 
strength of loss, bias and irregularity of the junction. The global attractor can 
contain co-existing local attractors e.g. a local chaotic attractor and a local 
regular attractor. In the infinite dimensional phase space setting, the bifur- 
cation is very complicated. Chaotic attractors can appear and disappear in a 
random fashion. Three types of attractors (chaos, breather, spatially uniform 
and temporally periodic attractor) are identified. The third topic is ratchet ef- 
fect. Ratchet effect can be achieved by a current bias field which corresponds 
to an asymmetric potential, in which case the flux dynamics is ever lasting 
chaotic. When the current bias field corresponds to a symmetric potential, 
the flux dynamics is often transiently chaotic, in which case the ratchet effect 
disappears after sufficiently long time. 
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1. Introduction 

The recently developed theory |12j on chaos in partial differential equations 
has the greatest potential of significance in its abundant applications in science and 
engineering. The variety of the specific problems also stimulates innovation of the 
theory. In these representative publications [13] [14] [15] [T7] [18] [19] [20] [2T]. 
two categories of the theory were developed. The category developed in [13] [14] 
|15j involves transversal homoclinic orbits, and a shadowing technique is used to 
prove the existence of chaos. This category is very complete. The category in [17] 
[18J [19| [20 [21 deals with Silnikov (non-transversal) homoclinic orbits, and a 
geometric construction of Smale horseshoes is employed. This category is not very 
complete. The main machineries for locating homoclinic orbits are (1). Darboux 
transformations, (2). Isospectral theory, (3). Persistence of invariant manifolds and 
Fcnichel fibers, (4). Melnikov analysis and shooting technique. Overall, the two 
categories of the theory on chaos in partial differential equations can be regarded as 
a new development at the intersection among Integrable Theory, Dynamical System 
Theory, and Partial Differential Equations [12] . In this article, we will apply the 
above chaos theory to study the chaotic flux dynamics in a long Josephson junction. 

A Josephson junction consists of three parts: two superconductors separated 
by a thin (~ 10 A) dielectric barrier. The main character of a Josephson junction 
is that it behaves as a single superconductor. In particular, electric flux called 
fluxons can travel through the dielectric barrier through Josephson tunneling. In 
general, the scenario is as follows: when two superconductors are separated by a 
macroscopic distance, their phases can change independently. As the two super- 
conductors are moved closer to about 30 A separation, quasiparticles can flow from 
one superconductor to the other by means of single electron tunneling. When the 
separation is reduced to 10 A, Cooper pairs can flow from one superconductor to the 
other (Josephson tunneling). In this case, phase correlation is realized between the 
two superconductors, and the whole Josephson junction behaves as a single super- 
conductor. This phenomenon was predicted by Brian Josephson in 1962 [9]. It has 
significant applications in quantum-mechanical circuits. The governing equations 
are [3] 

where V(t) and I(t) are the voltage and current across the Josephson junction, u 
is the phase difference between the wave functions in the two superconductors, the 
constant I c is called the critical current, and the constant ^ is the magnetic flux 
quantum. Interesting simple phenomena can be observed from the above governing 
equations, e.g. when V(t) is constant Vq in t, u is linear in t, then the current 
I it) will be an AC current with amplitude I c and frequency tVo- This shows that 
a Josephson junction can be a perfect voltage-to-frequency converter. Significant 
applications of a Josephson junction can be found in many areas, e.g. in medicine 
for measurement of small currents in the brain and the heart. Josephson junction 
may also provide key ingredients for future quantum computers. For more details 
on the physics of the Josephson junction, see [3]. 

When a Josephson junction has one or more dimensions longer than the rest, 
i.e. a long Josephson junction (LJJ), the phase difference u is also a function of 
spatial coordinates along the longer dimensions. In this article, we will focus on one 
longer dimension case. It is well known that the flux dynamics in a long Josephson 
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junction (LJJ) is described by the so-called perturbed sine-Gordon equation [3] 
(1.1) u tt = c 2 u xx + sinu + ef 

where again u is the phase difference between the two superconductors, and c is 
a constant. The above equation is the rescaled standard form. In dealing with 
real junctions, one must take into account losses, bias, and junction irregularities 
which influence the flux dynamics [3]. These effects are of perturbation nature, 
and accounted for by the term ef where e is the perturbation parameter. Different 
forms of ef can even by set up in experiments. Typically ef takes the following 
forms: 

• ef = —eau t + £7, where —eaut represents shunt loss and ej represents 
DC current bias [3]. 

• ef = —eau t + efiutxx + £7, where ej3ut X x represents longitudinal loss [3]. 

• ef = —ecxu t + £7 + easinwt, where ea sin cut represents AC current bias 

m- 

• ef — —eau t + £7 cos x + ea sin cot, where £7COSx represents AC current 
field [H [24]. 

• ef — —eaut + eg{x) + Easing, where eg(x) represents spatially periodic 
current field [4]. 

Of course, many other forms of ef can also be set up in experiments. There is an 
abundant literature on long Josephson junctions, for a sample, see [TJ [6J 1261 
l25ll23ll22lH0] . 



Travelling wave solutions of ( 1 . 1 1 satisfy an ordinary differential equation which 
had been studied numerically since as early as 1968 by Johnson [8] using an analog- 
digital computer, till most recently [5] via analytical and numerical tools. 

In this article, we shall study the full partial differential equation. When e = 0, 



equation (1.1 1 is the well-known sine-Gordon equation. Two well-known types of 
solutions to the sine-Gordon equation are the kink and breather solutions. The 
simplest kink solution is independent of the space variable x. Such a simple kink 
can be observed in experiemnts for Josephson junctions. As t — > ±00, the kink 
u(t) approaches e.g. and 2n. In terms of the current and voltage (I(t),V(t)) 
mentioned above, this corresponds to a closed loop — named a vortex or a fluxon. 
In the infinite dimensional phase space setting of the current article, this simple 
kink also represents a heteroclinic orbit. A breather oscillates periodically in time 
(breathing) and decays in space (x — > ±00). In our current setting, we switch the 
space and time, therefore, the breather is now spatially periodic and temporally ho- 
moclinic, i.e. a homoclinic orbit. For long Josephson junctions, spatial dependence 
is significant, dynamics beyond the simple kink becomes crucial. More sophisticated 
solutions incorporating both the kink and the breather have been constructed |16j . 
In the phase space, these represent heteroclinic orbits which depend on both time 
and space. Together, these heteroclinic orbits form a two-dimensional heteroclinic 
cycle. The current article will focus its study on the neighborhood of this hetero- 
clinic cycle. Three main topics will be covered. The first is chaotic flux dynamics 
and its prediction via Melnikov integrals (sections 4,5,6). The second is on the com- 
ponents of the global attractor and the bifurcation in the perturbation parameter 
(section 7). The third is ratchet effect (section 8). Section 9 is the conclusion. 
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2. The Phase Space 

We pose periodic boundary condition on the above perturbed sine-Gordon equa- 



tion (1.1 1 



(2.1) u(t,x + 2n) = u(t,x). 

This boundary condition corresponds to an annular LJJ |24j [4] , Convenient for 
mathematical studies, a more restricted condition — the even constraint may be 
imposed 

(2.2) u{t,-x)=u(t,x) 

which corresponds to a mirror symmetry of the annular LJJ. We restrict the pa- 
rameter c to the interval c £ (0.5,1) to minimize the number of unstable modes 



of u — to two. (1.1 1 is invariant under u —> u + 2tt. When e = 0, it is 
also invariant under u — > —u. ( |1.1| is globally well-posed |15j . i.e. for any 
(w°,it() G H n+1 x H n (the Sobolev spaces on [0,27r]) where n > 1 is an integer, 
there exists a unique mild solution to ( |l.l[ ), (u(t),u t (t)) S C([0, oo), x if"), 

such that (u(0),u t (0)) = ). One can introduce the evolution operator F* as 

(u(t),u t (t)) = For any fixed i G [0,oo), F* is a C°° map. The space 



H n+1 x if" will be the phase space where the perturbed sine-Gordon flow (1.1) is 
defined. 



3. Isospectral Integrable Theory 



When e = 0, the unstable manifolds of u + — and «_ = 2ir, W± form a 
2-dimensional heteroclinic cycle whose explicit expression can be obtained through 
Darboux transformations 



(3.1) 



(3.2) 



Ul 



U-2 



2arccos(tanhr) 

cr(± tanh t sech f cos x — sech r tanh f) 



-4 arctan 



1 — er(tanh r tanh f ± sech r sech f cos x) 
2tt — 2 arccos(tanh r) 

tr(± tanh r sech f cos a; + sech r tanh f ) 
1 — er(tanh r tanh f =p sech r sech f cos x) 



-4 arctan 



where the ranges of arccos and arctan are [0, 7r] and (— 7r/2, it/2), and (r, cr, f) are 
given by 



£ — p, cr = \fl 



at 



err — Ap, Ap = p — crp, 



and /J and p are two real parameters. This heteroclinic cycle satisfies both the 



periodic boundary condition and the even constraint \2.1\ and (2.2 1. A topological 
illustration of the heteroclinic cycle is shown in Figure |1] The heteroclinic cycle 
incorporates both kink and breather characteristics. The homoclinic degeneracy to 
be discussed later corresponds to the breather characteristic (when time and space 
are switched). The current article will focus its study near this heteroclinic cycle. 

Via isospectral theory, Melnikov vectors can be derived for the sine-Gordon 
equation |16j . Along each of u\ and u 2 in the above heteroclinic cycle, explicit ex- 
pressions of the Melnikov vectors can be obtained through Darboux transformations 
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(3.3) 
(3.4) 



dF 
du t 
dF 1 
du t 



647T ~ 

—r{G x 4> x + G 2 Af> 2 )(G 3 /0 1 + G 4 /0 2 ) , 



^(1-c 2 ) 
4c 



^l|$2| 2 )(^l|$l 



I/! | $ 2 1 



where Fq and i*i are two invariants of the sine-Gordon equation, vq — i/A, v\ 

-n/T^c 2 ], 
(3.5) = 2 

here = ir/2 for «i and 8 



cosh | cos | — i sinh | sin | 
— cosh | sin | — i sinh | cos | 



(3.6) 
and 

here £ = a; 



$ = 



7r/2 for 1*2; 



Vq4>\/4>2 



tp ■ 



cosh ^ cos I — i sinh ~ sin | 
— cosh | sin | — « sinh ^ cos | 



= ±7r/2, and the '±' corresponds to the '±' in u\ and 1*2; finally 



.^I^P + irrl^l 2 



|2 j_ ,,. <T» .1:2 
,2 7^2^ 



I/l|$l| 2 + I/l|$ 2 | 



G 2 

G 3 
G 4 



^ (^ 2 - t/i 2 )^>i$2 

PT|$i| 2 + j/i|$ 2 | 2 
z^o(^i -T7T 2 )¥7$ 2 
^il^iP + ^Tl^l 2 



= ^1 



| 2 ^l|$l| 2 + ^l|$ 2 

i/il$il 2 + ^r|$ 2 ; 



The two Melnikov vectors 



and will be used to heuristically measure the 



splitting of the heteroclinic cycle when the sine-Gordon equation is under pertur- 
bations ( 1.1 1. 
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3.1. Homoclinic Degeneracy . Inside the 2-dimensional unstable manifold 
0, there is a homoclinic orbit asymptotic to u + = 0. Let p — > — oo in 



Wl of u 



the expression of u% (3.1 ), we get 
(3.7) 



«1 



=F4 arctan 



a sech f cos x 



and U2 —* 2?r + U\. The other limit p 



1 — a tanh f ' 

> +oo produces equivalent expressions. 
Expression (3.7) represents a pair of homoclinic orbits asymptotic to u + = (cf: 
Figure 0. One can also generate the expression (|3.7|) directly from the Darboux 
transformation. 

In the limit of p — > — oo, the <f> defined in (3.5) has the limiting property: 



Thus the $ defined in (3.6 1 now has the form 



(3.8) 



i 

4 



_ 1 

4 



In the limit of p 



-oo, the Mclnikov vector §^ 

out 



du t 



(3.3), and the only Mel- 



nikov vector left is f^- which is given by (|3.4|) where <!> is now given by (|3.8|. This 



Melnikov vector is transversal to the pair of homoclinic orbits (3.7 1. From Figure 



[T] we see that the unstable and stable manifolds of u + = intersect along the pair 
of homoclinic orbits (3.7 1. Inside the d x = plane, these two manifolds intersect 
transversally at u + — 0, thus we can use to heuristically measure the split- 
ting of the homoclinic orbits when the sine-Gordon equation is under perturbations 



( 1 . 1 1 . The rigorous mathematical problems are often open 



4. The Case of / = —au t + f3u txx + 7 

4.1. Heteroclinic Cycle Search. In this case, we are interested in the possi- 
ble heteroclinic cycle for the perturbed sine-Gordon equation ( |1.1| , connecting the 
two fixed points 

u + = — arcsin(e7), u_ = 2tt — arcsin(e7). 

Such a possible heteroclinic cycle has the potential of inducing chaos. Unfortu- 
nately, as shown later, such a heteroclinic cycle does not exist. Nevertheless, we 
will present the failing arguments simply because the argument in the current case 
is the simplest. Arguments in other cases are often similar. Linearization at the 
fixed points leads to 

u tt = c 2 u xx + y/l - (ej) 2 u + e[-au t + (iu txx \. 

Let u = X)fc°=o u k(t) cos then 

u'i + e{a + f3k 2 )u' k + [c 2 k 2 - V 1 - (e7) 2 H = 0. 

Let Uk ~ e° fct , then 

fl k = ~e(a + 13k 2 ) ± \j ' y/l - (e 7 ) 2 - c 2 k 2 + Je 2 (a + (3k 2 ) 2 . 
The results of |20] |16j imply the following invariant manifold theorem. 
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Theorem 4.1. When c E (1/2,1) and e > is sufficiently small, each of 
the two fixed points u± has a 2- dimensional C m (m > 3) unstable manifold W± 
and a 2-codimensional C m center-stable manifold W± in the phase space (u, u t ) E 
jjn+i x jjn ( n > \) - When e > 0, each of u± has a 2-codimensional C m local stable 
manifold W± of size Q(\fe). In a 0(^/e) neighborhood ofu±, W± = W± . W± are 
C 1 in e for eE [0,e ) and some e > 0. When f3 > 0, at (u,u t ) E H n + 2 x H n + 2 , 
W± are C 1 in e for e € [0, e )- When [3 — 0, W± are always C 1 in e for e E [0, eo). 

To the leading order, the signed distances (which are certain coordinate differ- 
ences) between W± and W£? are given by the Melnikov integrals 

r+oo 



(4.1) 



M jt = 



27T (m 

1 du t 



-au t + I3u txx + 7] 



dxdt, 



where (J = 1,2; I — 0, 1) and w/s are given in (3.1 1- ( 3 . 2 ) . The common zero of the 
Melnikov integrals and implicit function theorem imply the intersection between 
W± and W^ s . The intersected orbits approach u± as t — > —00. As t — > +00, 
they can reach the 0(y/e) neighborhood of u=p where = W^ s , as shown in |20j . 
Thus, they approach as t — > +00, and form a heteroclinic cycle. So the key 
question now is whether or not the Melnikov integrals have a common zero. 

M jt = aM^ ] + [3M$ - 

where are functions of c and Ap = p — up, 

(dF e 



M 



(«) _ 



M 



C/3) 



M. 



it 
(7) 



+00 



+00 



2tt 



\ du t * 
(dF} 

1 o ^txx 

{ ou t 
du t J 



dxdt, 



dxdt, 



dxdt. 



It turns out from numerical calculations (we believe it is analytically provable) that 



all the M^f (z 

0,1) and M 1 ( J ) 
the system 



a, 13, 7; j = 1, 2; I = 0, 1) are real, and M 
= 0,1 



M 2 ( | } (z 



a,f3;£- 



Therefore, the common zero of Mj£ satisfies 



(4.2) 




f /3M^ 


+ 7 mW 


= 0, 


(4.3) 


aM[f - 


f 0Mff 


+ jM^ 


= 0, 


(4.4) 


«<' - 






= 0, 


(4.5) 


oM£> - 


f 


-1*41 


= 0. 



There is no non-trivial solution to this system (4.2 1-( 4.5 1 , which shows a failure in 
searching for a heteroclinic cycle. 

4.2. Heteroclinic Orbit Search . Next we try to search for an individual 
heteroclinic orbit. We only need to solve equations ( 4.2 1-( 4.3 1 or (4.4 1-( 4.5 ) . In 
cither case, we obtain the equation 



M {0) M h) 



s 
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Since both a and [3 have to be positive, x nas to be positive too. But direct 
calculation shows that x is always negative (see Table [TJ . This shows a failure in 
searching for even an individual heteroclinic orbit. 



Ap 


-10 


-5 





5 


10 


c = 0.55 


-1.0000 


-1.1617 


-1.4179 


-1.1617 


-1.0000 


c = 0.65 


-0.9998 


-1.0895 


-1.5753 


-1.0895 


-0.9998 


c = 0.75 


-0.9997 


-1.0423 


-1.7368 


-1.0423 


-0.9997 


c= 0.85 


-0.9990 


-1.0075 


-1.9272 


-1.0075 


-0.9990 


c = 0.95 


-0.9875 


-0.9870 


-2.0338 


-0.9870 


-0.9875 



Table 1. The table of x f° r the heteroclinic orbit case in section [4~2| 



4.3. Homoclinic Orbit Search. To search the homoclinic orbit asymptotic 
to u + — 0, we need to calculate the Melnikov integral: 

/ + OC rlTX t Qp 
J [d^H ^~ aUt + l3UtXX + ^ ^ 



where u\ is given in (3.7 1, t^- 

M = 

where 



is given by (|3.4| where $ is now given by (|3 . Sh 



Q 



M< a > = - 
AfW = 
AfW = 



+00 



2tt 



I du t 



Ui 



u I du t 
\du t 



dxdt, 



dxdt, 



dxdt. 



Numerical calculations show that M 1 - 1 ' is zero, and are of the same 

sign. So there is no nontrivial solution to M = 0. This shows the failure of search 
for a persistent homoclinic orbit. 

The failures of the Melnikov prediction on both a persistent heteroclinic orbit 
and a persistent homoclinic orbit indicate that there may be no chaos in the current 
case of perturbations. In general, it is difficult for autonomous perturbations as in 
the current case to generate chaos. That is, DC current bias usually does not 
generate chaotic flux dynamics. Our numerical simulations indicate that there is 
no chaos under the current perturbations. 

5. The Case of / = —au t + 7 + asinuit 
First consider the ordinary differential equation (ODE) by setting d x = 0, 
u tt — sin u + e [—au t + 7 + a sin ut] . 
The (e = 0) fixed point u = u t = turns into a periodic orbit when e > 0, 

= eui + e 2 u 2 + • • • . 
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To the leading order 0(e), 

Uitt = u\ + 7 + asirnvt 
which has the periodic solution 

a 

u\ = —7 smut. 

1 + LO 1 

Under the ODE flow, the periodic orbits and 27r+u» are normally hyperbolic and 
a Melnikov integral can locate a heteroclinic cycle connecting these two periodic 
orbits [16 . As a result, there is a region in the external parameter space that 
supports chaos |16j . 



Under the PDE flow ( 1.1 ), the stability of these periodic orbits and 2n+u* is 
complicated with parametric resonances as discussed in the open problems in |16j . 
As a result of this complication, the existence of a heteroclinic orbit connecting 
the two periodic orbits and the existence of chaos are open rigorous mathematical 
problems. Consider the extended system, 

u t t — c 2 u xx + sin u + e[— au t + 7 + a sin 9] , 



When e = 0, the periodic orbit u = Ut = 0, 9 6 T 1 has two unstable eigenvalues, 
two stable eigenvalues, and the rest neutral eigenvalues. This leads to the following 
invariant manifold theorem. 

Theorem 5.1. When c 6 (1/2, 1) and e > is sufficiently small, each of the 
two periodic orbits (u„,9 — Lot) and (2tt + u*,9 — cut) has a 2 co-dimensional 
C m (m > 3) center-unstable manifold W™, a 2 co-dimensional C m center-stable 
manifold W±, and a 4 co- dimensional C m center manifold W± in the phase space 
(u,Ut,6) G H n+1 x H n x T 1 (n>l). W£ = WJ" n W^ s . With W| as the base, 
W± and W± are fibered by 2-dimensional Fenichel fibers. 

Due to parametric resonances, dynamics inside W± still contains unknown num- 
ber of unstable and stable modes with growth or decay rates of 0(e). Melnikov 
integrals can detect orbits in W± that approach (2ir + it* , 9 = tot) or (u* , 9 = Lot) 
in backward time. In forward time, the destiny of these orbits is unknown — This 
is precisely the open problem |16j . 



5.1. Heteroclinic Cycle Type Connection Search. Specifically, the Mel- 
nikov integrals measuring the simultaneous intersection between W±" and W^f are 
given by 



M 



+00 p2n 



(9F lr . 

< — — -ati f + 7 + asmuit] 
I out 



dxdt 



r(s) 



10 



Z. C. FENG AND Y. CHARLES LI 



"0 



0, 1) and u/s are given in (3.1 1-(3.2|, are functions of c 



where (J = 1, 2; £ 

and Ap = p — <jp, M^f and also depend on to, and specifically 

»-foo />27r 



M 



(a) _ 



M 



(7) 



71/ 



M 



+ 00 



-OO 
+ 00 





2tt 



dxdr, 



o I <9u t J 

2tt 



dxdr, 



——sm{uT) 





cos(wr) 



dxdr, 



dxdr. 



It turns out from numerical calculations that Afj? and M$ are real and inde- 
pendent of j. The other integrals are the same as in previous section. Thus, the 
common zero of Mjg satisfies the system 



(5.1) 
(5.2) 
(5.3) 
(5.4) 



aM 



(a) 
10 



aM[f 



aM 
aM 



(a) 
10 
(a) 
11 



■ 7^' 

• 7 m&> 

7^' 



acos(u}p)M$ 
acos(ujp)M[f 
acos{u)p)M$ 
acos(ujp)M[f 



asm(u>p)M[^ > = 0, 
■aam(u)p)M$ = 0, 
■asm(up)M t $ = 0, 

asin^M^ = 0. 



There is no non-trivial solution to this system (5.1 1- (5.4 1, which shows a failure in 
searching for a simultaneous intersection between W±* and W^ s . Our numerical 
simulations indicate that there is no chaos of the type associated with heteroclinic 
cycles. 

5.2. Heteroclinic Orbit Type Connection Search . Next we shall search 
for an intersection between and WL S (or and W+) only, in which case 
we only need to solve equations ( |5.1[ )-(5.2 1. Eliminating a, we have 

j A + acos(Lop)B + asin(wp)C = 0, 

where 



A = 
B = 
C = 



(7) 

10 JW 11 
(a) 



M- 

M$M U 



(q) M (7) M (a) 
m u 1V1 W ' 



71/ 



M (c) M (a) 



, ' M y ■ 

'10 JW 11 

Thus as long as B and C are not simultaneously zero, there are always non-trivial 
solutions: 

a\/B 2 + C 2 ' 



sin(w/9 - 



where 



sin C/* = 



B 

So we get the criterion that when 



C 2 ' 



cosf* = 



C 

VB 2 + C 2 ' 



|ffl|>xl7l> where x = 



VB 2 + C 2 ' 
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Ap 


-10 


-5 





5 


10 


uj = 0.1 


1.01 


1.01 


0.95 


1.01 


1.01 


uj = 0.5 


1.32 


1.34 


0.58 


1.34 


1.32 


uj = 1.0 


2.51 


2.53 


0.67 


2.53 


2.51 


uj = 2.0 


11.61 


11.68 


10.25 


11.68 


11.61 


uj = 3.0 


55.74 


56.00 


53.06 


56.00 


55.74 


uj = 4.0 


264.3 


262.0 


149.7 


262.0 


264.3 


uj = 5.0 


1290. 


1231. 


507. 


1231. 


1290. 



Table 2. The table of \ f° r the hctcroclinic orbit case in subsec- 
tion [O] c= 0.75. 



there is an intersection between Wf" and Wt s (or and W+). The values of x 
is shown in Table [2] 

We conducted a numerical simulation of the dynamics with the following setup: 

e = 0.1, a = 1.0, uj = 0.55, c = 0.75, a = 0.257, 7 = 0.5; 

and 64 elements are used to divide the spatial period [0,27r], the time step is 1/40 



of the forcing period finally the initial condition is given by (3.1 1 with (t = 0, 
p = 0, p = 0, hence Ap = 0). The numerical result indicates chaotic dynamics 
as shown in Figures |2][3] One can see a clear monotonely shifting in the u(x, t) 
plot. This indicates that the chaos is due to the persistent heteroclinic orbit (the 
hetroclinic cycle does not persist as shown in last subsection). The shifting is due 
to the fact that when it ends up outside the "eye" in Figure [TJ the orbit will travel 
up to the other eye. If we view the range of u(x,t) as a circle (i.e. mod 27r), then 
the shifting will disappear. The u t plot in Figure [3] also illustrates this effect. 



Figure 2. The chaotic dynamics in subsection 



5.2 



(it- plot). 
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Figure 3. The chaotic dynamics in subsection 5.2 (u t -plot). 



5.3. Homoclinic Type Connection Search. To search for the homoclinic 
type connection near u + = 0, we need to calculate the Melnikov integral (which 
measures the intersection between and WSf): 

■-oo />2?r 



(5.5) 



M 



oo JO 



dF 1 
du t 



-au t + 7 + a sin tot] 



dxdt, 



where u\ is given in (3.7 1, is given by (3.4 1 where $ is now given by (3.8). 



M = aM {a) + 7 M (7) + acos(ujp/(T)M (c) + asin(wp/cr)M (s 



where 



Af<"> 



+00 /"27T 



-oo JO 

OO />27T 



-oo ^0 

+ 00 p27T 



dxdf, 



oo Jo 
-f oo p2ir 



{dF\ 
du t 



I 5iT C0S ^ T / cr ) 



dxdr, 



dxdr, 



dxdf. 



As before, = 0. It turns out that both and are zero (in fact 

has zero spatial mean). Thus there is no nontrivial solution to M = 0. 
This shows a failure in searching for a persistent homoclinic orbit. Our numerical 
simulations indicate that there is no chaos of the type associated with homoclinic 
orbits. 
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6. The Case of / = —au t + 7 cos a; + asinwi 

In this case, there is no guarantee that the (e = 0) fixed point u = u t = 
will persist into a periodic orbit. The existence of heteroclinic orbits and chaos is 
a far more open problem than the case in the previous section. Again consider the 
extended system, 



u tt = c 2 u xx + sin it + e[—au t + 7 cos a; + asinO], 

= LU. 



When e = 0, the periodic orbit u — u t — 0, 9 e T 1 has two unstable eigenvalues, 
two stable eigenvalues, and the rest neutral eigenvalues. This leads to the following 
invariant manifold theorem. 

Theorem 6.1. When c € (1/2,1) and e > is sufficiently small, there are a 
2 co-dimensional C m (m > 3) center-unstable manifold W± , a 2 co- dimensional 
C m center-stable manifold W±, and a 4 co- dimensional C m center manifold W± 
in the neighborhoods ofu= 0,27r; ut = 0, 6 G T 1 (where u = corresponds to '+' 
and u = 2ir corresponds to '—') in the phase space (u,u t ,0) G H n+l x H n x T 1 



(n>l). W. 



± 



n Wl s . With Wl as the base, Wl v and W^ s are fibered by 



2-dimensional Fenich el fibers. 



6.1. Heteroclinic Cycle Type Connection Search . Mclnikov integrals 
can detect orbits in n W^ s . Unfortunately both the forward and the backward 
destinies of such orbits are unknown. The Melnikov integrals are given by 



Mi 



-aut + 7 cos x + a sin cot] 



' +oc f 2 * (dF) 
'-co ■/,. I , „-„, 
= aM ( .f + jM^ + a cos(ojp)M ( j f + a sin(wp)M ( ; ) , 



dxdt 



where (j = 1,2; i = 0, 1) and u/s are given in (3.1 1-(3.2 1, M$ are functions of c 



and Ap — p — ap, M^f and M^f also depend on to, and specifically 



M {a) 



M 



(7) 



j i' 



M, 



M 



v <- 27T (dF e 

I du t Ut 
+00 r**,Qfi 



— OO 

OO /■27T 



1 9u t 
dF, 



00 jo { 9u t 

+ OO /-27T I Qp f 

00 Jo 



\ du t 



cos 2 
sui(u;t) 
cos(wr) 



dxdr, 



dxdr, 



dxdr, 



dxdr. 



It turns out from numerical calculations that all the integrals Af -f , M 



r(7) 



A/ 



(s) 

and -M-^ are real and independent of j. The other integrals are the same as in 
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previous section. Thus, the common zero of Mj£ satisfies the system 



(6.1) 
(6.2) 
(6.3) 
(6.4) 



aM 



(a) 
10 



aM[f 



1< 7) 

■ 7 m&> 
• 7 m&> 



acos(wp)M 1 ^ 
acos(wp)M 1 ^ 
acos(wp)M 1 ^ 
acos(wp)M 1 ^- ) 



asm^^-M^ = 0, 
•osiu^Mff = 0, 
asin(wp)M-[Q' 1 = 0, 
asin(a;p)M}^ = 0. 



Notice that equations ( 6.3 )- ( 6.4 1 are identical with equations (6.1 1-( 6.2 ) . This de 



generacy prohibits the application of the implicit function theorem in measuring 
the intersection between W± u and . So we have no conclusion as whether or 
not and W^f can simultaneously intersect. 

We conducted a numerical simulation of the dynamics with the following setup: 



3.0; 



e = 0.1, a = 5.1, u = 1.0, c = 0.75, a = 0.257, 7 

and 64 elements are used to divide the spatial period [0, 2tt], the ti me st ep is 1/40 of 
the forcing period ^ — 2tt, finally the initial condition is given by (3.1 1 with (t = 0, 
p = 0, p = 0, hence Ap = 0). The numerical result indicates chaotic dynamics as 
shown in Figure [4] In fact, it seems that the chaotic dynamics loops around the 
eye in Figure [l] i.e. around a cycle, rather than monotonely travels up to different 
eyes as in Figure [2] Possible explanation is that the chaotic dynamics is induced 
by persistent heteroclinic cycles even though our Melnikov calculation above has 
no conclusion. 




Figure 4. The chaotic dynamics in subsection |6.1 



6.2. Heteroclinic Orbit Type Connection Search . Next we shall search 
for an intersection between and WL S (or Wf" and W+) only, in which case 
we only need to solve equations ( |6.1[ )-(6.2 1. Eliminating a, we have 

j A + acos{ujp)B + asin(wp)C = 0, 
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Ap 


-10 


-5 





5 


10 


LO = 0.1 


36.05 


1.27 


2.35 


1.27 


36.05 


lo = 0.5 


47.17 


1.67 


1.44 


1.67 


47.17 


LO = 1.0 


89.36 


3.16 


1.64 


3.16 


89.36 


lo = 2.0 


413.3 


14.59 


25.23 


14.59 


413.3 


lo = 3.0 


1985. 


69.93 


105.9 


69.93 


1985. 


lo = 4.0 


9411. 


327.2 


387.8 


327.2 


9411. 


lo = 5.0 


45950. 


1538 


1214. 


1538. 


45950. 



Table 3. The table of x for the heteroclinic orbit case in subsec- 



wherc 



A = 






B = 




-m£>m$ 


C = 







Thus as long as B and C are not simultaneously zero, there are always non-trivial 
solutions: 

jA 

sm(Lop + 6i) = 



where 



VB 2 + C 2 ' 



B C 
s'mdi = =, cos 0i = 



VB 2 + C 2 ' yjB 2 + C 2 

So we get the criterion that when 

I A I 

(6.5) \a\ > xh\, where \ 



VWTC 2 ' 



there is a heteroclinic orbit. The values of x 1S shown in T able [3] The chaotic 
dynamics in Figure E] is in the parameter regime predicted by (6.5 1 



6.3. Homoclinic Type Connection Search . To search for the homoclinic 
type connection near u + = 0, we need to calculate the Melnikov integral (which 
measures the intersection between and 



(6.6) 



M 



+00 p2-7T 



du t 



-aut + 7 cos x + a sin Lot] 



dxdt. 



where u\ is given in (3.7 1, is given by (3.4 1 where $ is now given by (3.8). 



M = aM {a) + 7 M (7) + acos(Lop/(j)M (c) + asin{ujp/ a)h&\ 
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where 



M (a) 
M (c) 



■ r2v f dF x 
\ du t 

I du t 



iii 



oo JO 

OO />27T 



oo JO 

-f-OO /•2-7T 



dxdr, 



dxdf, 



oo Jo 

+ 00 /•27T 



00 JO 



9^i . , , 
- — sm(u;T/crJ 
out 

|_ cosW(7) 



dxdf, 



dxdr. 



As before, and M( 5 ) are zero, but now M^ 1 ' is not zero. So Af = implies 

that 

7 = X a, where x = - M (a) / M {l) . 
When c = 0.75, we found that 



X 



-6762.7. 



That is, comparing to the shunt loss amplitude a, the AC current field amplitude 
7 needs to be quite large to generate possible homoclinic chaotic current (CC). We 



Figure 5. The chaotic dynamics in subsection |6.3| 

conducted a numerical simulation of the dynamics with the following setup: 

e = 0.01, a = 0, uj = 0.55, c = 0.75, a = 0.001, 7 = -6.7627; 

and 32 elements are used to divide the spatial period [0,27r], the time step is 1/80 
of the forcing period ^j, finally the initial condition is given by (3.1 1 with (t = 0, 
p = 0, p = 0, hence Ap = 0). The numerical result indicates chaotic dynamics 
as shown in Figure [5] This chaos is still jumping around the heteroclinic loops, 
nevertheless most of the time it loops very close to the homoclinic orbits. 
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7. The Components of the Global Attractor and the Bifurcation in the 

Perturbation Parameter 



Here we study the case of / = — aut + 7 + a smut in (1.1 1 which was also 
studied in section [5j The setup is as follows: 

(7.1) = 1.0, ui = 0.55,c=0.75, a = 0.257, 7 =0.5; 

and 64 elements are used to divide the spatial period [0,27r], the time step is 1/40 
of the forcing period — . 



The perturbed sine-Gordon system (1.1 1 has a global attractor up to the trans- 
lation u — > u + 2ir [7] . Here we are interested in the detailed structure of the global 
attractor. In particular, we are interested in the local attractors inside the global 
attractor. Components other than the local attractors, often can only attract mea- 
sure zero sets of initial conditions. Figure [6] shows two local attractors inside the 
global attractor when e = 0.1. For the initial condition (it*, u* t ) given by (3.1 1 with 
t = 0, p = 0, p = (hence Ap = 0, and it* = 77); the solution approaches the 
chaotic attractor Figure [6j a). For the initial condition 

(tl, tit) = (77,0.1 + 0.01r(ar)) 

where r{x) is a random function depending on x\ the solution approaches the pe- 
riodic attractor Figure Hb) . To explore the initial conditions more systematically, 




(a) 



(b) 



Figure 6. Two local attractors inside the global attractor. (a) is 
a chaotic attractor, and (b) is a periodic attractor. 

we use the following homotopic initial conditions 

(it, lit) = (1 - A)(u*,i4) + A{u, u t ) 

where A G [0, 1] and (it, Ut) = (77,0.1) i.e. we dropped the random perturbation in 
(it, iii), which does not seems to affect the final attractor. Another topic that we are 
interested in is the bifurcation in the perturbation parameter e. By combining the 
two parameters e and A, we obtain the component-bifurcation diagram (Table [4]), 
where H is the heteroclinic orbit given by (3.1 1, Q is the quasiperiodic solution of 
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e 


A = 0.00 


A = 0.05 


A = 0.10 


A = 0.15 


A = 0.20 





H 


Q 


Q 


Q 


Q 


0.01 


U 




u 


u 


u 


0.02 


U 


U 


u 


u 


u 


0.03 


B+ 


u 


u 


u 


u 


0.04 


C 


c 


c 


u 


u 


0.05 


c 


c 


u 


u 


u 


0.06 


u 


c 


u 


u 


u 


0.07 


c 


c 


u 


u 


u 


0.08 


c 


u 


c 


u 


u 


0.09 


u 


c 


c 


u 


u 


0.1 


c 


u 


u 


u 


u 



Table 4. The component-bifurcation diagram in section [7] H: 
heteroclinic orbit, Q: quasiperiodic orbit, U: spatially uniform 
and temporally periodic attractor, £>+,£>_: breather attractors, 
C: chaotic attractor. 



the integrable sine-Gordon equation (e — 0). U is spatially uniform and temporally 
periodic attractor as depicted in Figure [6jb) . It is originated from the integrable 
steady state (u, u t ) = (it, 0) and modulated by the temporally periodic forcing 
asinut. B + is a breather attractor as depicted in Figure [7J It is originated from 
the center of a loop of one of the small figure-eights in Figure [l] and modulated by 
the perturbation. B_ is a breather attractor similar to B + except that the hump 
is located at the spatial periodic boundary. It is orginated from the center of the 
other loop of the small figure-eight. C is the chaotic attractor as depicted in Figure 
|6ja). When the value of e or A is larger than those in Table [4] the attractor is 
always U. From Table [4] one can see that there is no clear boundary between chaos 
attractor and regular attractor in the (e, A)-planc. The chaotic attractors appear 
and disappear in an irregular fashion when the parameters are varying. In the 
perturbation parameter e, the bifurcation does not follow any simple bifurcation 
paradigm in low dimensional systems. 

8. Ratchet Effects 

The long Josephson vortex ratchet effect refers to the effect of nonzero temporal 
average of the voltage output Ut when the temporal part of the input forcing / 



(1.1 1 has a zero temporal average. There has been a lot of recent interest on the 



Josephson vortex ratchet, see [4] and the references therein. Here we will study the 



following form of the input forcing / ( 1 . 1 1 



(8.1) / = — au t + g(x) + asin ut, 

where g(x) is a spatially periodic current field of zero mean. It turns out that the 
existence of the ratchet effect depends on whether or not the spatial potential G 
(G x = g) is symmetric as observed experimentally [4]. 
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Figure 7. The breather attractor B + in Table [I] 



First we study the asymmetric potential case: 



g(x) = 7 for \x — n\ < 



7 



for other x £ [0, 2ir]. 



48' 47 

Figure [8] corresponds to the following setup: 

(8.2) e = 0.1, a = 2.0, w = 0.55, c = 0.75, a = 0.14, 7 = 100; 

and 64 elements are used to divide the spatial period [0, 2tt], the time step is 1/40 of 
the forcing period finally the initial condition is given by (3.1 ) with (t = 0, p = 0, 
p = 0, hence A/? = 0); the temporal average of u t is done over 256 data points, i.e. 
6.4 [256/40] times of the forcing period — . One can see a clear nonzero temporal 
average of the voltage Ut at x = 0. That is, there is a clear ratchet effect. The 
spatio-temporal profile of u t is clearly chaotic in time, while the spatio-temporal 
profile of u is also chaotic and dominated by drifting in time. 

Next we study the symmetric potential case: g(x) — 7 cos a; which is the case 
studied in Section [6] Figure [9] corresponds to the same setup as (8.2 1 and the rest. 



One can see that the temporal average of the voltage Ut at x — is approaching zero 
in time. That is, there is no ratchet effect in long term. The spatio-temporal profile 
of u t is transiently chaotic in time, while the spatio-temporal profile of u is also 
transiently chaotic in time and its drifting in time is mild. Figure [10] corresponds 
to the same setup as (8.2 1 and the rest except a = 5.0. One can see that the 

is also slowly approaching zero in time. 
The spatio-temporal profile of Ut is 



temporal average of the voltage u± at x 
That is, there is no ratchet effect in long term, 
transiently chaotic in time with a longer term, while the spatio-temporal profile of 
u is also transiently chaotic in time and it almost has no drifting in time. 
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(c) 

Figure 8. The temporal average of Ut at x = is shown in (a), 
while (b) and (c) are spatial-temporal profiles. This is the ratchet 
effect in the asymmetric potential case. 

9. Conclusion and Discussion 

Via a combination of numerical and Melnikov integral studies, we find that 
DC current bias cannot induce chaotic flux dynamics, while AC current bias can. 
The existence of a common root to the Melnikov integrals is a necessary condition 
for the existence of chaotic flux dynamics. The global attractor can contain co- 
existing local attractors e.g. a local chaotic attractor and a local regular attractor. 
In the infinite dimensional phase space setting, the bifurcation is very complicated. 
Chaotic attractors can appear and disappear in a random fashion. In the parameter 
space, there is no clear regular boundary between local chaotic attractors and local 
regular attractors. Three types of attractors (chaos, breather, spatially uniform 
and temporally periodic attractor) are identified. Ratchet effect can be achieved 
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(c) 

Figure 9. The temporal average of Ut at x = is shown in (a), 
while (b) and (c) are spatial-temporal profiles. This corresponds 
to no ratchet effect in the symmetric potential case. 



by a current bias field which corresponds to an asymmetric potential as observed in 
experiments [4], in which case the flux dynamics is ever lasting chaotic. When the 
current bias field corresponds to a symmetric potential, the flux dynamics is often 
transiently chaotic, in which case the ratchet effect disappears after sufficiently long 
time. 

Due to its infinite dimensionality, numerically exploring the entire phase space 
is impossible. Here we focus upon an interesting neighborhood. It is entirely 
possible that other novel structures are hidden somewhere else in the phase space. 
Also due to the infinite dimensionality, the link between existence of chaos (and 
homoclinic orbit) and existence of a common root to the Melnikov integrals becomes 
weaker. Finally, chaos in the infinite diemsional phase space is often transient. After 
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(c) 

Figure 10. The temporal average of Ut at x = is shown in (a), 
while (b) and (c) are spatial-temporal profiles. This also corre- 
sponds to no ratchet effect in the symmetric potential case. The 
difference from Figure [9] is that the value of the parameter a is 
different. 



sufficiently long time, the seemingly chaotic dynamics may converge to a regular 
attractor. 

Inside the global attractor, there may be many invariant components. An in- 
teresting topic is to classify these components. Of particular interest are those 
components which are local attractors. Some local attractors may be chaotic, while 
others may be regular. Different initial conditions may lead to different local at- 
tractors. A complete classification of all these local attractors is very challenging, 
especially in the infinite dimensional setting. This article only explores a selected 
set of initial conditions. When the value of the perturbation parameter changes, the 
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dynamics undergoes bifurcations. Unlike in low dimensional systems, bifurcations 
in the infinite dimensional system are much more complicated. 
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